Designing Robust Unitary Gates: Application to Concatenated Composite Pulse 
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We propose a simple formalism to design unitary gates robust against given systematic errors. 
This formalism generalizes our previous observation [Y. Kondo and M. Bando, J. Phys. Soc. Jpn. 
80, 054002 (2011)] that vanishing dynamical phase in some composite gates is essential to suppress 
pulse-length errors. By employing our formalism, we derive a new composite unitary gate which 
can be seen as a concatenation of two known composite unitary operations. The obtained unitary 
gate has high fidelity over a wider range of error strengths compared to existing composite gates. 
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I. INTRODUCTION 

Noise and errors are obstacles against reliable control 
of a quantum system. Noise, i.e., random unwanted dis- 
turbance to a quantum system we concern, has been at- 
tracting much attention of theoreticians. Many ideas to 
suppress noise have been proposed [l|, in quantum com- 
putation, which requires precise control of quantum sys- 
tems 0. Geometric quantum gates (GQGs) [34l2j|. that 
are based on holonomy [l3l421j . are such examples. On 
the other hand, errors, i.e. systematic imperfection in 
control parameters, have also been attracting attention, 
due to their importance in realistic situations. 

To tackle the latter problem, one may decompose a 
given unitary gate into a sequence of several unitary oper- 
ations, whose time-ordered product reproduces the given 
unitary gate [23 - l32j . Then the sequence becomes robust 
against given systematic errors by tuning the parameters 
in the constituent unitary operations. Such sequences for 
a two-level system is well-known as composite pulses in 
NMR and have been designed by employing various tech- 
niques, such as the Magnus expansion [25j j and quater- 
nion algebra [30| . for example. In the following, we often 
use a "composite gate" to denote a composite pulse when 
it is regarded as a quantum gate. 

As mentioned in [23[ , there are several lines of thought 
to understand composite pulses in a unified manner. Mo- 
tivated by these, we have proved in Q that GQGs for a 
two-level system are insensitive to an error in the am- 
plitude of the control parameters and shown that many 
existing composite pulses are regarded as GQGs. This 
shows that we can coherently interpret several composite 
pulses for a specific systematic error in terms of geometric 
phases. In this paper, we extend our former work [9] in 
order to include general systematic errors. The derived 
conditions are simple enough to be understood straight- 
forwardly and applicable not only to GQGs, but also to 
gates involving dynamical phases. Our new formalism is 
applicable straightforwardly to multi-partite systems ei- 
ther. As a demonstration of our formalism, we design a 
new composite pulse robust against the most important 
systematic errors in NMR. The obtained pulse sequence 



can be seen as a concatenation of two composite pulses 
derived in [30| and has high fidelity over a wide range in 
the error parameter space. This pulse sequence cannot 
be constructed by iterative expansion [24], H(| . 

This paper is organized as follows. In Sec. II, we 
present our formalism to design unitary gates robust 
against general errors. The robustness of the GQGs 
against the pulse-length error is generalized to arbi- 
trary non-degenerate multi-level systems in the contin- 
uous time cyclic evolution, which results in Abelian geo- 
metric phases. Further, we derive robustness conditions 
systematically based on our theory, after which discrete 
time formalism is introduced. In Sec. Ill, the developed 
formalism is applied to construct concatenated pulse se- 
quences, which are robust against the most important 
systematic errors in a two-level system. Sec. IV is de- 
voted to conclusion and discussions. 



II. ROBUSTNESS CONDITION 

First, we define robustness of a gate. Consider the 
special unitary group SU(n), whose dimension as a group 
manifold is N := n 2 — 1, and its Lie algebra Su(n). Wc 
introduce a complete set of orthogonal Hermitian basis 
(generators) of su(n), {r^ | r M = rLfi — 1, . . . , Y} with 
respect to the Hilbert-Schmidt inner product, and a time- 
dependent real Y-vector X(t) = (Ai(i), . . . , Ajv(i)) m the 
Y-dimensional parameter manifold Ai. By choosing a 
continuous path X(t) in yVf, we define a family of time- 
dependent Hamiltonians 



H(\(tj) := A^i)^ 6 su(n), 
whose time-evolution operator 



U x (t,0) :=Texp 



daff(A(*)) 



(1) 



(2) 



is an element of SU(n). Here, T denotes the time-ordered 
product and we set h — 1 and employed the Einstein 
summation convention for Greek indices. U\{t, 0) reduces 
to the identity operator 1 at t = 0. Note that the system 
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is a two-level system when n = 2 and t p = er M /2, where 
cr„ is the //-th component of the Pauli matrices. 

Now, let us scale t G [0, 1] and require that the time 
evolution operator U\(1,0) implements a target gate U 
at t = 1; C/a(1,0) = £/. We define J7a(1,0) is a robust 
gate if the condition 



U x+sx (l,0) = U + O(\5X\ 2 ) 



(3) 



is satisfied for a given SX(t) with |<5A(t)| -C |A(f)| for ev- 
ery t. To find the robustness condition, let us rewrite 
the LHS of Eq. ([3]) in the interaction picture. Consider 
the dynamics under a Hamiltonian H(X(t) + SX(t)) = 
H(X(t)) + H(5X(tj), where H(5X(t)) is regarded as a per- 
turbation. Then, we obtain 



C/a+5a(1,0) = C/C// a (1,0), 



where 



[// A (l,0) = Te X p 



dtH!(SX(t)) 



(4) 



(5) 



is defined through the interaction picture Hamiltonian 
Hj(6X(t)) := U x (t,0)iH(6X(t))U x (t,0). (6) 

Equation Q requires that U\(1,0) is robust against a 
given error 5X(t) if J7/ A (l,0) = t. 

To proceed further, let us rewrite Ug X (l, 0) in the RHS 
of Eq. Q as the Dyson series. Then, we obtain 



U x+Sx (l, 0) = U- iUAW + 0(|5A J | 2 ), 



where 



AW 



dtff/(5A(t)) 



(7) 



(8) 



and SX 1 is defined implicitly as a solution of Hi(SX) = 
H(SX r ). Since 0(|<5AT) = C>(|<5A| 2 ) from Eq. @, 
C/ A (l, 0) is robust against given 8X to 0(|<5A| 2 ) if and only 
if 



AW = 0. 



(9) 



Note that the robustness condition (J9)) is derived with- 
out assuming an explicit form of the Hamiltonian and is 
applicable to many physical systems. 



Classification 



So far, we have considered general 5X(t). From now 
on, we restrict ourselves within systematic (determinis- 
tic) errors for SX(t). By definition of a systematic error, 
SX(t) takes a form 



SX p (t) = F„(A(f)), 



(10) 



where F — . . . , F^) is an unknown vector function 

defined in M. We assume the RHS of Eq. (fl"Q|) admits 
an expansion 

F»(X(t)) = + f^X v (t) + f^ P X v {t)X p (t) + ■■■ ,(11) 



where / M , f pu , fpup, ■ ■ ■ are constant tensors. Hereafter 
we do not write the time-dependence of functions explic- 
itly to simplify equations. Substitution of Eq. (fTTj) into 
Eq. (JHJ) shows that Eq. © is satisfied for any F if 

dtTp = 0, / dtTpX„ = 0, / dtfpX v X p = 0, . . .(12) 



(i 



hold simultaneously. Here, we utilized the expression (TTJ) 
for H(SX(t)) and introduced generators fp(t) in the in- 
teraction picture: 



,(t) := £/ A (i,0)S?7 A (i,0). 



(13) 



The first condition in Eq. (|12p requires that the effect of 
a constant error F^ = / M on the time-evolution operator 
vanishes. 

We turn to the second condition in Eq. (TT2")) . To find its 
implication, f pu is decomposed into the sum of tensors, 
each of which is an irreducible representation of su(n) : 

fttv = S^fpp/N + (fpu - /^)/ 2 

+ [(U„ + U)/2 - 5p V f pp /N] . (14) 

Employing an analogy to fluid mechanics, each irre- 
ducible tensor in the RHS can be thought of as an error 
which causes a uniform expansion, rotation and torsion of 
X(t), respectively. According to the decomposition (|14p . 
we find 



dtfp,Xp= / dtU x (t ) 0^H(X(t))U x {t,0) = 0,(15a) 
Jo 

[ dt (fp,X„ — T v Xp) = 0, (15b) 
Jo 

dt [(f M A„ + f v Xp) /2 - S^fpX p /N} = (15c) 

as the sufficient conditions for the robustness against the 
corresponding errors. 



B. Geometric Phase Gate and Norm Error 
Compensation 

Let us consider the case when there is only a norm 
error in A exists. In this case, Eq. (jlip is reduced to 



Fp.(X(t)) = f PP Xp(t), 



(16) 



and we have to consider only f dtf p X p in Eq. (|15a|) for 
evaluating its robustness. 

Let us introduce \tp a (t)) ■= U\(t, 0)\tf) a (0)) and as- 
sume non-degeneracy of eigenvalues throughout time- 
evolution. By taking the expectation value of Eq. (|15a[) 
with respect to \ip a (0)}, we find that the dynamical 
phases [151 ] 



7d 



dt(Mt)\H(Ht))\Mt)) 



(17) 
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must vanish for all a. Now let us consider a case in which 
V'a(O)) is a cyclic state of U\(1,0) with a phase 7° e R 
(l5j . that is, IV'a(O)) is an eigenvector of £/a(1, 0) with the 
eigenvalue e 17 ° , 

|Va(l)> = (7 A (1,0)|V^(0)) =c i7 > a (0)). (18) 

We have the spectral decomposition of U\(l, 0) in terms 
of these mutually orthogonal {|*0o(O))}a=l as 

n 

[/ A (l,0)=^e^|^ a (0))(^ a (0)|. (19) 

a=l 

Let us recall that a cyclic state admits the Aharonov- 
Anandan phase [l5[ 

7g Q := 7 a -7d- (20) 

Thus, we realize a non-trivial unitary gate (7° ^ 0, 
mod 2ir for some a) robust against the error on the norm 
of the vector X(t), if it has a nonvanishing geometric con- 
tribution 7° under the condition (I15a|l . which leads to 
7j = for all a. This observation confirms that our for- 
malism is a proper continuous time generalization of the 
previous work [9j, which revealed that composite pulses 
robust against the pulse-length error are the GQGs. 

C. Discretization 

Next, let us divide the temporal interval [0, 1] into k 
intervals, in each of which the Hamiltonian is constant. 
With this piecewise constant Hamiltonian, time-ordering 
in the time-evolution operator U\(t, 0) is simply an or- 
dered product of time-evolution operators defined for 
each time-independent Hamiltonian. This means that we 
restrict ourselves within the gate U which is decomposed 
into a product 

U = U xk (t k ,t k - 1 )---U x i(t 1 ,t ), (21) 

where = to < t\ < 4 • • < tfc-i < tfc = 1 and X 1 is the set 
of constant parameters in the Hamiltonian H(X(ti)) cor- 
responding to the interval [ij_i,tj]. Then, by definition, 
the unitary operator in the presence of an error 5X(t) has 
a similar decomposition 

U\+s\(l, 0) = U X k + s\k (t k , tk-i) ■ ■ ■ U\i + $\i (t 1, t ), (22) 

where S X 1 is the deviation of the parameters in the inter- 
val [£{_!,£{]. Note that 8X l is independent of time when 
A* itself is independent of time. We write U\i(ti,ti~i) as 

c/*(<i,ti_i) -> : = cx p {-' im W) > ( 23 ) 

where m l := X l (ti — tj_i). By denoting the correspond- 
ing error as dm 1 , the unitary operator Ux*+s\* (*i> is 
replaced similarly as 

Uxi+sx* (ti, U-i) -v ^(^ 4 + «m*). (24) 

Let us introduce 

1=772^ and JW* := ^m^, (25) 



and recall the well-known formula 

e A+B = [ dx5(l - x)e Ax 
Jo 

+ [ dx [ dyS(l -x- y)e Ax Be Ay + ■■■ (26) 
Jo Jo 

for matrices A and B of the same dimension. Then, by 
setting A = —iW l and B — — iSW 1 and neglecting higher 
order terms with respect to B, we obtain 

R(m l + 6m l ) » R(m l ) - ii?(m l )<5PF;, (27) 

where 

6Wj:= [ dxe lxW *5W l e-' lxWl . (28) 

Thus, it follows that AW defined in Eq. ||7J) is given by 

k 

AW = ^V i - 1 hw i I V i - x , (29) 
»=i 

with 

=ii(m*)---i?(m 1 ) for i = 1,2, . . ., (30) 

and V° = 1. Note that we utilized the identity 
R(m k ) ■ ■ ■ i?(m i ) = UV 1 - 1 " 1 to derive Eq. Char- 
acterization of error (14) is applied under this discretiza- 
tion. 

In case [W 1 , SW l ] = for all i, we easily find 

5m% = em;, (31) 

which is the error generating the uniform expansion of 
the norm of the parameter vector m M . Then, we observe 
further simplification of AW: 

k 

Aff = e^y M V¥- 1 . (32) 
i=l 

This is nothing but the central quantity considered in Q , 
in which the usefulness of Eq. (|3"2"|) has been presented in 
detail. 



III. APPLICATION: CONCATENATED 
COMPOSITE PULSE 

As an example of our discretization formalism, with 
NMR and similar systems in mind, we construct an SU(2) 
gate robust against two errors defined by the following 
error operator 

5W l = tW l + e'\m l \a z /2. (33) 

The first term in the RHS causes the pulse-length error 
while the second term causes the off-resonance error in 
the terminology adopted from NMR. The error operator 
(|33p simultaneously represents the most important errors 
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inherent in quantum control in NMR and a system with 
an analogous Hamiltonian. 

Bearing the situation in NMR setups in mind, where 
there is no ct 3 in the Hamiltonian ([1} , we construct gates 
robust against errors (|33l) under the restriction 



m l = 0j(cos</>j,sin<^j,O), 



9i > 0. 



(34) 



Now, we concatenate two composite gates, each of 
which is composed of three simple unitary operations 
of the form ([23")) . One is the pulse sequence called 
Compensation for Off-Resonance with a Pulse SEquence 
(CORPSE), which is robust against the off- resonance er- 
ror and consists of three commutative pulses (30j . In 
other words, given a target U = R(m) with m = 
0(cos(j), sincj), 0), the elementary pulses in CORPSE are 
given by 



(35) 



and 



91 =6</2-K + 2ni7r, 

9 2 = -2k + 2n 2 7r, 

9 3 = 9/2 -k + 2n 3 n, 



where rii, ri2, n.3 G Z are chosen so that 6$ > and 

sin(0/2) 



re := arcsm ■ 



(36) 



(37) 



This pulse was obtained in [30( by making use of quater- 
nion algebra. Its robustness against random telegraph 
noise is examined in [3ll ] in comparison with other pulse 
sequences. CORPSE is used to compensate for the off- 
resonance error in NMR experiment [29(. 

The other composite pulse is the SCROFULOUS, 
which is an acronym of Short Composite Rotation For 
Undoing Length Over and Under Shoot. SCROFULOUS 
is robust against the pulse-length error and made of a n- 
pulse sandwiched between two identical pulses |30( • More 
precisely, three elementary pulses satisfy 



0i cos ( 



2 )+vr/2 = O. (38) 



SCROFULOUS is a generalization of a composite gate 
proposed in [33[, which has experimental confirmation 
for its robustness [29| . 

Now we are ready to construct the concatenated com- 
posite pulse by combining three CORPSE gates to form 
a SCROFULOUS gate. First, from Eq. (251), we find 



5Wj = eW l + e'Rim^as sin(0</2). (39) 



for the error operator (|33|) . For this, Eq. (|29|) for k = 3 
gives 

AW = R{m 1 )' { R{m 2 ^5WfR{m 2 )R{m 1 ) 

+ Rimy SWj Rim 1 ) + 8W} . (40) 

For notational convenience, we work with 

UAW = ei?(m 3 )S + e'T, (41) 



instead of AW, where 

S = W 3 R(m 2 )R{m 1 ) + Rim 2 )W 2 Rim 1 ) 

+ Rim^Rim^W 1 (42) 

and 

T = sm(e 3 /2)a 3 R(m?)R(rn l ) + sm{e 2 /2)R{m 3 )u 3 R{m 1 ) 

+ sin(6»i/2)i?(m 3 )i?(m 2 )CT 3 . (43) 

For CORPSE, we observe that T — as expected, since 
it is designed so as to compensate for the off-resonance 
error. Further, from Eq. (|3"5|) . we have \W l , R(mP)\ = 0, 
which leads to 

AW = eiW 1 + W 2 + W 3 ) (44) 
= e [9/2 + (n-i - n 2 + n 3 )n] (cos^cti + sin</>cr 2 ). 

Let us choose {n-i} so that they satisfy 

ni - n 2 + 713 = 0. (45) 

Then we have 



AW = em,,T,, 



(46) 



which is nothing but the pulse-length error acting on the 
target unitary operator. This clearly tells us that we can 
compensate for both systematic errors simultaneously, if 
we use a Concatenated Composite Pulse (CCCP) se- 
quence of three CORPSE sequences under the condition 
that they compose the SCROFULOUS when combined 
together. We call this concatenated pulse by CORPSE 
In SCROFULOUS-CCCP, or CIS-CCCP for short, in the 
following. 

One could alternatively try a concatenation of three 
SCROFULOUS pulses under the condition that they 
compose the CORPSE. This pulse sequence is, however, 
not robust in the sense of Eq. ([3]): Each constituent 
SCROFULOUS in the pulse sequence leads to S = but 
T / a 3 sin(0/2). This implies AW ± e / E/tcr 3 sin(0/2), 
that is, in view of the second term in the RHS of Eq. (|39p . 
the error which is not compensated for by each SCRO- 
FULOUS pulse cannot be regarded as the off-resonance 
error and the overall CORPSE fails to eliminate it. 

We would like to emphasize that the CIS-CCCP se- 
quence cannot be generated by iterative expansions 
[24l |26| . An iterative expansion is composed of consecu- 
tive applications of various pulse sequences, each of which 
is created from a given pulse sequence by i) a permuta- 
tion and ii) a shift of the rotation axes in the xy-plane of 
constituent pulses. Then, one cannot create a CORPSE 
with the total rotation angle 7r by operations i) and ii) 
on a generic CORPSE. This proves impossibility of de- 
signing the CIS-CCCP by iterative expansions. 

It is of interest to compare the fidelity of CIS-CCCP 
sequence with those of CORPSE and SCROFULOUS. 
The fidelity with respect to the target unitary gate U 
is defined by the absolute value of the Hilbert-Schmidt 
inner product: 



1 



Tr [U^R{m k + 5m k ) ■ ■ ■ Rim 1 + 8m 1 )] . (47) 
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(a) Plain (b) CORPSE (c) SCROFULOUS (d) Concatenated pulses 

FIG. 1. Fidelity T of (a) a plain pulse (b) CORPSE (m = n 3 = 1, n 2 = 2), (c) SCROFULOUS and (d) the concatenated 
pulse sequence as a function of the error strength constants e for the pulse-length error and e for the off-resonance error. The 
target unitary operator U is exp(— ma2/2). In whiter area the fidelity T yields higher value. 



Note that k = 3 for CORPSE and SCROFULOUS, 
whereas k = 9 for CIS-CCCP. 

Our interest lies in the weak error strengths region 
—0.1 < e,e' < 0.1, since the accuracy threshold theorem 
requires the error probability less than O(10 -3 ) for fault 
tolerant quantum computation [3j. From Fig. 1, we im- 
mediately observe two features of the CIS-CCCP. First, 
the CIS-CCCP has characteristics of both the SCRO- 
FULOUS and the CORPSE pulses, as expected: The 
CIS-CCCP is robust along the lines ee' = 0, whereas the 
CORPSE is robust along e = and the SCROFULOUS 
along e' = 0. Second, the whiter area, the higher fidelity 
region, of the CIS-CCCP is considerably wider than those 
of the CORPSE and the SCROFULOUS combined to- 
gether. This observation indicates that the concatena- 
tion of composite pulses results in an even more robust 
pulse sequence. 

In closing this section, let us show the difference be- 
tween the CIS-CCCP and the composite pulse proposed 
by Alway and Jones (32[. Their composite pulse also 
compensates the pulse-length error and the off-resonance 
error simultaneously, but implements only 7r-pulses along 
an axis in the ccy-plane (see Appendix [AJ. Here, we 
should note that any combinations of 7r-pulses along axes 
on the xy-plane are reduced to either 7r-pulse on the same 
plane or the pulse along the z-axis: 



R(m k ) ■ ■ ■ R(m x ) 



R(m), 
R(m z ), 



1,3,. 
2,4,. 



(48) 



Here, m? = 7r(cos<^ 7) sin0j,O) and m = 7r(cos <f>, sin^, 0), 
where <\> is read as a function of {<j)j}- We introduced 

to 2 = (0, 0, 20) with = Ej=i(02j - faj-i + tt). Equa- 
tion (|48l) can be derived by mathematical induction with 
respect to k. Therefore, it is impossible to implement ar- 
bitrary one-qubit rotations as combinations of their com- 
posite pulse sequences. This clearly shows that we cannot 
realize universal gate set by using their composite pulse 
sequence. In contrast, the CIS-CCCP does not have such 
a restriction as Eq. (|48p. and implements any one-qubit 
unitary operation by using Euler angles. 



IV. CONCLUSION AND DISCUSSIONS 

In this paper, we proposed a simple formalism to de- 
sign unitary gates robust against systematic errors whose 
magnitude are unknown. By using this method, we sys- 
tematically derived various criteria which admitted lucid 
interpretations. 

We designed a pulse sequence robust against two types 
of systematic errors (|3"3"]) simultaneously. We design a 
new composite pulse that is the SCROFULOUS out of 
the CORPSEs in order to take advantages of these com- 
posite pulses. Our approach is straightforward than the 
quaternion algebra that relies on brute force calculation. 
Our pulse sequence has controllable free parameters; cf. 
the pulse sequence proposed in (32|, which is also robust 
against the errors (|33[) , but implements only 7r-pulse gate 
on the xy-plane. In contrast, our pulse sequence real- 
izes arbitrary one-qubit unitary gate robust against com- 
bined errors (|33[) and will find an important application 
in implementation of universal gate set out of low qual- 
ity gates. These features show the usefulness of our for- 
malism as a guiding principle to construct unitary gates 
robust under coexisting deterministic errors. 

We would like to stress that the our scheme is applica- 
ble not only to NMR, but also to other physical systems 
since the condition ^ is formulated independently of the 
Hamiltonian. For example, the effective Hamiltonian of a 
quantronium superconducting qubit takes the NMR form 

m 



H = —hvRo(cos^ai + sin £er 2 ) 



(49) 



with an additional off- resonance error term hAv<j3, 
where h is a constant and Av is the detuning. It was 
demonstrated in [34{ that CORPSE indeed suppresses 
the off-resonance error. A similar demonstration of the 
effectiveness of CORPSE has been made for a neutral 
atom qubit to suppress effective microwave detuning 
across the qubit ensemble [Hj]. Our concatenated com- 
posite pulse is applicable to these systems. 

Our method is also useful for designing a robust two- 
qubit gate. Nonetheless, it requires intensive analytical 
as well as numerical analysis and is beyond the scope of 
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the present paper. Our preliminary result shows that a 
two-qubit gate with error in the coupling strength be- 
tween qubits may be made robust against the error by 
decomposing the gate into a relatively small number of 
gates, which will be reported elsewhere. 
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Appendix A: Alway and Jones' gate 

In Sec. IIIIl we discussed the composite Tr-pulse de- 
signed by Alway and Jones [32J , which compensates both 
the pulse-length error and the off-resonance error. Sup- 
pose we are to implement a pulse sequence whose rota- 
tion axis is in the xy-pl&ne. Then the pulse sequence 
employed in [32j is robust under the error (j33|) if and 
only if the rotation angle is tt. 

The "if" part is proved by Alway and Jones [32[ as 
mentioned above. Let us briefly reproduce their result 
to establish notation and conventions. For notational 
convenience, we introduce U(0) := exp (— \0o\/2). (We 
consider arbitrary rotation axes at the end of Appendix.) 
Given the target U(tt), we can design the following pulse 
sequence: 

U oi U p iU(tt), (Al) 
where we introduced two partial sequences 

3 7 

f/pi = n ^k)' u o* = n i? ( m ')' ( A2 ) 

i=l i=4 

whose elementary pulses are given as 
9 2 = 2tt, 9i = ir (i?2) 

01 = 03 — <f>ph 02 = 30 p l, 

4> A = 2lT - 6 = 7T - or , 5 = -07 = -0 o r, (A3) 



with p i = or = arccosj — 1/4). U P \U(%) is a composite 
pulse known as BB1 [3g, which is robust against the 
pulse-length error, while U or U(ir) is a composite pulse 
robust against the off-resonance error. Note that 



U p i = U OI = t 



(A4) 



in the absence of errors. This is the sequence proposed 
in j32|. 

Now let us prove the "only if" part. Suppose we want 
to implement U(9), which is robust against simultaneous 
errors by employing the sequence (|A1|) . Here we assume 
9 is not fixed to tt and p i and or are adjusted so as to 
make the gate robust for a given 8. After simple calcula- 
tion, we obtain the zeroth and the first order error terms 
as 



U{e)-ie-a 1 U{e)-ie'sm-a 3 , 



(A5) 



for the target U{6). In the same way, we have 

[/pi = 1 — 27riecos0pi(Ji (A6) 

and 

U OI = 1 + 4ie' cos or o"2- (A7) 

Taking product of (|A6]l and (TAT)) and evaluating 

the coefficients of e and e', we find that this composite 
pulse is robust against the simultaneous errors only if 9, 
p i and or satisfy the conditions: 

9 9 
27rcos0 pl +- = 0, 4 cos or cr 2 £/■(#)— sm-a- 3 = 0. (A8) 

These conditions have a unique nontrivial solution 

P i = 0or = arccos(— 1/4), 9 = it, (A9) 

which shows that the rotation angle of the target pulse 
must be tt. 

Clearly the rotation axis can be any direction in the 
xy-plane by a simple redefinition of the coordinate axes 
or by applying the similarity transformations around the 
z-axis to each constituent pulses. This completes the 
proof of the "only if" part. 
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